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Abstract 



We compute lattice renormalisation constants of local bilinear quark operators 
for overlap fermions and improved gauge actions. Among the actions we consider 
are the Symanzik, Liischer-Weisz, Iwasaki and DBW2 gauge actions. The results 
are given for a variety of p parameters. We show how to apply mean field (tadpole) 
improvement to overlap fermions. The question, what is a good gauge action, is 
discussed from the perturbative point of view. Finally, we show analytically that 
the gauge dependent part of the self-energy and the amputated Green functions are 
independent of the lattice fermion representation, using either Wilson or overlap 
fermions. 
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1 Introduction 



Lattice calculations at small quark masses require an action with good chiral properties. 
The same is true for calculations of matrix elements of certain operators, which otherwise 
mix with operators of opposite chirality. Ginsparg- Wilson fermions [Tj have an exact chiral 
symmetry on the lattice [2J, and thus are well suited for these tasks. A further advantage 
is that they are automatically 0(a) improved 0- Overlap fermions |U E] provide a 
four- dimensional realisation of Ginsparg- Wilson fermions. The massless Neuberger-Dirac 
operator is defined by 



where Dw is the Wilson-Dirac operator, and p is a real parameter corresponding to a 
negative mass term. At tree level p must lie in the range < p < 2r, r being the 
Wilson parameter. Exact chiral symmetry on the lattice has its price however. Numerical 
implementations of the overlap operator are significantly more expensive than Wilson 
fermions. In spite of this difficulty, calculations with overlap fermions are progressing 
rapidly, and we expect to see many more results in the near future. 

It is important to use a good gauge field action. The cost of the overlap operator is 
governed largely by the condition number of X'X. This number is greatly reduced for 
improved gauge field actions [7j. For example, for the tadpole improved Luscher-Weisz 
action we found a reduction factor of > 3 compared to the Wilson gauge field action jH]. 
The reason is that the Luscher-Weisz action, as many other improved actions j^j, sup- 
presses unphysical zero modes, sometimes called dislocations. A reduction in the number 
of small modes of X^X appears also to result in an improvement of the locality of the 
overlap operator [7]. 

To obtain continuum results from lattice calculations of hadron matrix elements, the 
underlying operators have to be renormalised. A perturbative calculation of the corre- 
sponding renormalisation constants is always the first step. Recently, the renormalisation 
constants of bilinear [TDJ HU E2| and four-quark [T3] operators have been computed to 
one-loop order for overlap fermions and Wilson gauge field action. Our aim is to extend 
the calculation to general, improved gauge field actions with up to six links, including the 
tree-level Symanzik action, the Liischer-Weisz action, the Iwasaki action, and the DBW2 
action. In this paper we shall consider local bilinear quark operators first. Operators in- 
cluding covariant derivatives will be treated in a separate publication. Preliminary results 
of the calculation have been reported in [S]. 

The paper is organised as follows. In Section 2 and Appendix A we give the details of 
the calculation. The basic results are given in Section 3. In Appendix B we show that the 
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gauge dependent terms in the Green functions are not only independent of the gauge field 
action, but also of the type of fermion. Tadpole improvement is an important issue in 
lattice perturbation theory. In Section 4 and Appendix C we tadpole improve our results. 
Finally, in Section 5 we give our conclusions. 



2 Calculational details 



We denote the lattice action by 

S = So + Sf , (2) 

where Sq is the gauge field action, and Sp is the fermion action. Let us discuss the gauge 
field action first. We consider the following class of improved actions: 
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where &pi aq uette is the standard plaquette, while the remaining C/'s cover all possible closed 
loops of link matrices of length six along the edges of the hypercubes, as indicated in Fig.QJ 
It is customary to impose the normalisation condition 



c + 8ci + 16c 2 + 8c 3 = 1 



(4) 



In Section 0] we will also consider actions which fulfill this condition only in the limit 
g — > 0, which is sufficient to ensure the correct continuum limit. Calling the coefficient in 
front of the plaquette part the lattice coupling /3, we have the obvious relation 
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If the improvement is performed only on-shell, one of the six-link contributions can be 
set to zero ^3]- In general, the parameters q are chosen according to an approximate 
renormalisation group (RG) analysis or using tadpole improved perturbation theory. 



For perturbation theory expansions, we write the SU(3) link matrices as 
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(6) 



where = T a A'J 1 , T a being the generators of the Lie algebra, and perform a weak 
coupling expansion in g. The calculations are carried out in a general covariant gauge, 
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Figure 1: The meaning of the four terms in the gauge field action J3J). 



distinguished by the gauge parameter £ = 1 — a. Landau gauge corresponds to £ = 1, 
while the Feynman gauge corresponds to £ = 0. The gluon propagator (in the infinite 
volume) is obtained from the lowest order expansion of (jHj): 
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±sin^, P = J2K- 0) 



The coefficients C$ are related to the coefficients q of the improved action by 

Co = c + 8ci + 16c 2 + 8c 3 , Ci = c 2 + c 3 , C 2 = ci - c 2 - c 3 . (10) 
(We will usually work with actions in which Cq = 1; see eq. (JU).) 
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In lattice momentum space the gluon propagator D )iv {k) is given by the set of linear 
equations 

2_j G^ p (k) — - — - kjzp 



p 



D pv {k) = 8^ . (11) 



For the plaquette action the gluon propagator reads 




ku k v 



r;;;; ; —(/,) -(^-^] ■ (12) 

Expressions for the more general action can be found in the literature ^HJ QHj- We 
will use dimensional regularisation to regularise the loop integrals. Thus we need to 
know Dp V for arbitrary dimensions. In general, can be given in closed form only for 
integer dimensions. Only for the special case C2 = can we derive explicit expressions 
for arbitrary dimensions. A way out is to split the gluon propagator into two parts: a 
singular part, which can easily be extended to arbitrary dimensions, and a finite part, 
which does not need to be regularised. This is achieved by writing 

Dp U = C l D^ ette + A£V . (13) 

It is readily seen that Dp V and Cq 1 D^ uette have the same infrared singularity. The 
finite part of the gluon propagator AD^ is obtained by solving (fTTj) for the difference 
Dfj, v — Cq 1 Dv^"i nette (in four dimensions). It can be written 

4 4 4-n 

azv = v Yl D -( k > c ) K n + E E D ^ c ) k l m+1 k l n+l ■ ( 14 ) 

n=0 n=0 m=0 

The scalars D n and D mn = D nm are rational functions of k^ and Cj. They are listed in 
Appendix A for the case C = 1. 

The final results cannot be expressed in analytic form (as a function of C$) anymore. 
We therefore have to make a choice concerning the parameters. We restrict ourselves to 
the most popular actions in this class: 

Tree-level Symanzik 

Cl = -1/12, C2 = c 3 = 0. (15) 
Tadpole improved Liischer-Weisz (TILW) [H Q3J EH] 

c 2 = 0. (16) 
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Once we have chosen /3, the other parameters are fixed [T§] : 
(1 + 0.4805 a) on 0.03325 a 
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3.06839 



We have picked the following values 
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Co 
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(17) 



(19) 





Cl 


c 3 


8.60 


-0.151791 


-0.0128098 


8.45 


-0.154846 


-0.0134070 


8.30 


-0.159128 


-0.0142442 


8.20 


-0.161827 


-0.0147710 


8.10 


-0.165353 


-0.0154645 


8.00 


-0.169805 


-0.0163414 



which corresponds to lattice spacings a = 0.084 — 0.136 fm. These j3 values are currently 
being employed in quenched Monte Carlo simulations, or most probably will be in the 
near future. 

Iwasaki 1 



Ci = —0.331, c 2 = c 3 = 



(20) 



DBW2 



ci = -1.4086, c 2 = c 3 
In all these cases c is chosen to satisfy eq. 



0. 



(21) 



Which action is the best one? They all reduce or eliminate 0(a 2 ) corrections - at the 
tree level, the one-loop level, or beyond. In Section 4 we compare the actions with regard 
to their perturbative merits. 
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The action for massless overlap fermions is given by 

S F = TpD N ip . (22) 

The Neuberger-Dirac operator was already given in eq. (0). The Wilson-Dirac oper- 
ator Dw reads 

~ k(V* + V„) - arV^V J , (23) 



D 



w 

where V M is the lattice forward covariant derivative: 

V^O) = - \u x ^i)(x + afi) -il>(x)] . (24) 
a L J 

For < p < 2r the correct spectrum of massless fermions without doubling is obtained [0] . 

The lattice Feynman rules for overlap fermions, being originally derived in [221 EH], are 
collected in Appendix A. The regularisation of the infrared divergences follows [21], which 
was adapted from [22]- The calculations are done analytically as far as possible. To do so, 
we have extended our Mathematics programme package, which we developed originally 
for Wilson [21] and clover fermions [26] . to overlap fermions with improved gauge actions. 

We performed several tests to check the code. Operators with gamma structures 
T and must give the same result. We have performed calculations independently 
for both sets of gamma matrices. Although the analytic expressions looked very much 
different, we found complete agreement. Another condition is that all terms oc 1/a cancel. 
Furthermore, it is expected [TT] that the gauge dependent parts of the one-loop lattice 
integrals are independent of the particular choice of fermion propagator (be it Wilson or 
overlap), as a consequence of the Ward identities. This has been verified numerically and 
analytically. In Appendix B we give an analytic proof. 



3 Renormalisation 



To obtain finite answers, the lattice operators 0(a) must, in general, be renormalised. 
Ignoring operator mixing, we define renormalised operators by 

O s (p) = Z*(a^)0(a), (25) 

where S denotes the renormalisation scheme. The renormalisation constants Zq are often 
defined in the MOM scheme first by computing the gauge fixed quark propagator Sn 
and the amputated Green function A© of the operator O: 

Zf OM {a^)S N \ 2=n2 = 5 tree , (26) 



A a 2= 2 — + other Dirac structures . (27) 





p p + k p p p 

Figure 2: One- loop diagrams contributing to the quark self-energy. 



(Note that = 1/Z 2 .) The renormalisation constants can be converted to the MS 
scheme, 

Zli(a, l*) = zS MOM Z*?™(a, n) , (28) 



,MOM ■ lit 



where Z^f%' M M is calculable in continuum perturbation theory. 



We shall restrict all our numerical calculations to the case r = 1. The optimal choice 
for p appears to be p ~ 1.4. We will present results for p = 1.3, 1.35, 1.4, 1.45 and 1.5, 
which should cover the most interesting region. Any other value in this region may be 
obtained by inter- or extrapolation. 

Self-energy and wave function renormalisation 

Let us consider the massless quark propagator Sn first. The inverse of Sn can be 
written 

^ X = iy(l- ^fZi) , (29) 

with Cf = 4/3. The diagrams that contribute to Ei to one-loop order are shown in Fig. 
The integral to be evaluated is |21| 

Sl / -j^Y^^P + k) (3°) 



16vr 2 



Explicit expressions for the propagator and vertex functions are given in Appendix A. 
Putting everything together, we finally obtain 

Ei(a,p) = (1-0 log(aV)+4.7920U + &s. (31) 

The results for are given in Table ^ for the gauge field actions and p parameters 
mentioned. The result for the plaquette action (cq = 1, C\ — c 2 — C3 = 0) agrees with 
Refs. [El [TT] . The gauge dependent part of (|3~Tj) is the same for all actions: 

(F - lE + !)£ = 4.79201 £, (32) 



8 



Action p 







1 QD 
l.oU 


l.oo 


i a n 
JL.4U 


1 A K 
1.40 


1 KCt 
JL.OU 


Plaquette 




-26.400 


-25.031 


-23.766 


-22.592 


-21.501 


Symanzik 




-21.774 


-20.713 


-19.732 


-18.823 


-17.978 


TILW, /3=t 


3.60 


-19.216 


-18.329 


-17.509 


-16.749 


-16.044 


TILW, f3=i 


3.45 


-19.124 


-18.243 


-17.429 


-16.675 


-15.974 


TILW, f3=i 


3.30 


-18.997 


-18.124 


-17.319 


-16.572 


-15.878 


TILW, p=t 


3.20 


-18.918 


-18.051 


-17.25 


-16.508 


-15.819 


TILW, p=l 


5.10 


-18.817 


-17.957 


-17.163 


-16.426 


-15.743 


TILW, p=l 


5.00 


-18.692 


-17.841 


-17.054 


-16.325 


-15.648 


Iwasaki 




-15.960 


-15.295 


-14.681 


-14.112 


-13.584 


DBW2 




-10.373 


-10.096 


-9.841 


-9.605 


-9.386 



Table 1: The contribution 6s to the self-energy for the various actions and parameters, 
with 

F = 4.3692252338748 , -f E = 0.57721566 . (33) 

We explain the reason for this observation in Appendix B. For the wave function renor- 
malisation constant we then obtain 



rMOM ( „ \ _ i 9 



2/T 

1 - f-4 [2(1 - logM + 4.79201 i + 6s] 

167T Z 



(34) 



In the MS scheme this becomes 



2 X~Y 

Zf^(a, /i) = 1 - f-4 [2(1 - logM + 3-79201 £ + 6 S + 1] • (35) 



For Wilson fermions (with r = 1), on the other hand, we have 

^SiisonK AO = 1 - TTTT [2(1 - Mo//) + 3.79201 £ + 12.8524] . (36) 



We observe that the constant, gauge independent terms in eqs. fjHHjl and (J3HJ1 have opposite 
sign. So for \i = 1/a the overlap Z$ is greater than one, while Z</,,wiison is less than one. 
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Figure 3: One-loop vertex diagram contributing to the amputated Green functions. 

Quark bilinears 

Let us consider local operators of the form 

O x = ^(x)T x ^(x) (37) 

with 

r s = i, r p = 75 , r v = 7 „, r A = 7 M 75, r T = a^ 75 , (38) 

i.e. X = S, P, V, A and T. As the operators are local, operator tadpole and cockscomb 
diagrams do not contribute [21]. This leaves us to compute the vertex diagram shown in 
Fig. El We denote the amputated Green function of the operator Ox by A x . Thus we 
have to evaluate the integral 



(39) 



x Y x S N (p + k) V lv (p + k, k) D^(k) . 

The final result is 

A s > p = {1, 75} + ^£ [(£ - 4) log(oV) " 5-79201^ + b s , P ] {1, 75} , (40) 

1071" 

= {7M,7M7b} + ^{7M[(e-l)log(«V)-4.7920ie + M 

-2(1-0^/) {1,75}, (41) 
Aj„ = <V7b + [f log(aV) " 3.79201^ + frr] <v 7 fi . (42) 



10 



Action p 







l.oll 


1 Q K 
l.OO 


i a n 
1.4U 


1 A K 
1.40 


1 KC\ 
1.0U 


Plaquette 




11.083 


11.214 


11.343 


11.472 


11.600 


Symanzik 




10.621 


10.739 


10.856 


10.972 


11.086 


r 1 ITT TT T ,0 

TILW, 13= 


8.60 


10.275 


10.385 


10.493 


10.599 


10.705 


TILW, 0= 


8.45 


10.261 


10.370 


10.478 


10.584 


10.689 


TILW, 0= 


8.30 


10.241 


10.350 


10.457 


10.563 


10.667 


TILW, 0= 


8.20 


10.229 


10.337 


10.444 


10.549 


10.653 


TILW, /3= 


8.10 


10.213 


10.321 


10.427 


10.532 


10.636 


TILW, p= 


8.00 


10.192 


10.300 


10.406 


10.510 


10.613 


Iwasaki 




9.641 


9.736 


9.830 


9.922 


10.012 


DBW2 




7.457 


7.515 


7.571 


7.626 


7.679 



Table 2: The contribution bs,p to A 5 ' p for various actions and parameters. 



Action p 







1.30 


1.35 


1.40 


1.45 


1.50 


Plaquette 




6.343 


6.345 


6.348 


6.351 


6.354 


Symanzik 




6.332 


6.333 


6.335 


6.338 


6.340 


TILW, p= 


8.60 


6.325 


6.327 


6.329 


6.330 


6.332 


TILW, p= 


8.45 


6.325 


6.327 


6.328 


6.330 


6.332 


TILW, p= 


8.30 


6.325 


6.326 


6.328 


6.330 


6.332 


TILW, p= 


8.20 


6.325 


6.326 


6.328 


6.330 


6.331 


TILW, p= 


8.10 


6.324 


6.326 


6.328 


6.329 


6.331 


TILW, p= 


8.00 


6.324 


6.326 


6.327 


6.329 


6.331 


Iwasaki 




6.316 


6.318 


6.319 


6.320 


6.321 


DBW2 




6.302 


6.302 


6.303 


6.303 


6.304 



Table 3: The contribution by,A to A V ' A for various actions and parameters. 
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Action p 







1.30 


1.35 


1.40 


1.45 


1.50 


Plaquette 




4.096 


4.056 


4.016 


3.977 


3.938 


Symanzik 




3.958 


3.914 


3.871 


3.828 


3.785 


TILW, 0=f 


3.60 


3.851 


3.804 


3.759 


3.713 


3.668 


TILW, /3=S 


3.45 


3.846 


3.800 


3.754 


3.708 


3.664 


TILW, f3=i 


3.30 


3.840 


3.793 


3.747 


3.702 


3.657 


TILW, /3=* 


3.20 


3.836 


3.789 


3.743 


3.698 


3.652 


TILW, /3=S 


3.10 


3.831 


3.784 


3.738 


3.692 


3.647 


TILW, (3=i 


3.00 


3.825 


3.778 


3.731 


3.686 


3.640 


Iwasaki 




3.651 


3.601 


3.551 


3.501 


3.452 


DBW2 




2.943 


2.881 


2.819 


2.759 


2.698 



Table 4: The contribution fry to A T for various actions and parameters. 
The finite terms bs,p, b v< A and &t are collected in Tables |2l El and EJ For the renor- 



malisation constants we then obtain, using (|27jl and ([34)1 . 

Zf° M {a,n) = l-^[-61ogM-^& s ,p + M, (43) 

^TW) = l-f^(V,A + & S ), (44) 

Zf OM (a,/i) = l-^[21og(a/i)+£ + & r + 6 s ] . (45) 

In the MS scheme the renormalisation constants read 

Zf*{a, /i) = 1 - f—J [-6 logM - 5 + b s , P + 6 E ] , (46) 

Z™(a,ii) = 1-^(V,a + &e), (47) 

Z^(a )A i) = 1 - [2 log(a/x) + 1 + &r + to] . (48) 



The conversion factors Z ' are universal, and are the same as for the plaquette 
action and Wilson fermions [21]. 
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4 Tadpole improvement 



The appearance of gluon tadpoles in lattice perturbation theory make the bare coupling 
constant g into a poor expansion parameter. It was proposed [2Z| that the perturbative 
series should be rearranged in order to get rid of these contributions. Tadpole improvement 
is a technique for summing, to all orders, the numerically large perturbative contributions 
arising from tadpole diagrams. It is implemented by a mean field renormalisation of 
the link matrices, U x ^ — *> U XiI1 /uq, where uq is the mean value of the link, defined to 
be the fourth root of the expectation value of the plaquette 'measured' in Monte Carlo 
simulations, as given in eq. I)18j) . In this Section we apply tadpole improvement, better 
called mean field improvement, to our results, and to operators involving overlap fermions 
in general. 

Our situation is more complicated than the standard case for two reasons: firstly we 
are using overlap fermions rather than Wilson or clover fermions, and secondly we are 
using gauge actions which are more complicated than the basic plaquette action. Both 
these changes require some thought. 

First, let us run briefly through the procedure for the case of Wilson fermions. Tadpole 
improved renormalisation constants are defined by 



where Zq is the mean field approximation of Zq, while the second factor on the r.h.s. 
is computed in perturbation theory. We find Zef Wilson by looking at the operator Green 
function and quark propagator in the mean field approximation. The mean field result 
for the amputated Green function for an operator with rip covariant derivatives is 



This does not depend on the choice of fermion action. In the mean field approximation 
the inverse quark propagator for massless Wilson fermions is 




(49) 





(50) 




(51) 




(52) 




(53) 
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The result changes when we consider overlap fermions. The mean field result for the 
amputated Green functions, eq. (joTIj) . is unchanged. However, for overlap fermions is 
no longer given by (|52j) . which will lead to changes in Zq as well. To find Z^ we need to 
calculate the overlap operator in the mean field approximation. The mean field result for 
the Wilson-Dirac operator in momentum space is 



U W 



- I u Q ^ 1 7^ sin a Pv + 4r - u r ^ cos a Pn ) 

= i $ Uq H (1 - u ) + -ap 2 u + 0(a 2 ) . 

a 2 

Inserting this result into the Neuberger-Dirac operator (^J) gives 



(54) 



nMF _ P U 

N 



p - 4r(l - u ) 
We can now compute Z¥ F : 



+ -ap — r + 0(a 

2 p-4r(l-« ) 



(55) 



r p — 4r(l — Mo) 
If we combine this result with eq. (j5U|) . we finally obtain 

C = . (57) 

° p-4r(l-u ) 1 ; 

It is required that p > 4r(l — u ). The reason for this inequality is easy to understand. 
Reverting temporarily to writing the fermion matrix in terms of a hopping parameter k, 
the p definition in eq. (0) is equivalent to p = 4r — 1/ (2k). To really have a negative mass 
in X, so that the overlap procedure works, requires k > n c , which implies p > 4r — 1/(2k c ). 
In the mean field approximation k c = l/(8rw ), leading to the inequality p > 4r(l — u ). 
In other words, the additive renormaliation of mass, which occurs with Wilson fermions, 
means that in the interacting case p > is not enough to cause X to have a negative 
mass. We need p to be large enough to overcome this additional mass. 

Note that Z^f F will be larger than 1 for reasonable values of p, while in the Wilson 
case Z ¥yyuson < 1- O ur mean field calculation gives an explanation of the observation 
in Section 3, that the one- loop b^s have the opposite sign for overlap fermions than for 
Wilson fermions. 

Let us now turn to the calculation of the second, perturbative factor on the r.h.s. 
of eq. (}4T?|) . Having removed the large tadpole contributions, we need to re-express the 
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Action 


a 
P 


h 


K u 


u 


Plaquette 


6.0 


7T 2 


7T 2 


0.59368 


Symanzik 




0.732524 7T 2 






TILW 


8.60 


0.590078 7T 2 


0.549643 vr 2 


0.66018 


TILW 


8.45 


0.585039 7T 2 


0.543338 vr 2 


0.65176 


TILW 


8.30 


0.578127 7T 2 


0.534971 vr 2 


0.64252 


TILW 


8.20 


0.573860 7T 2 


0.529705 vr 2 


0.63599 


TILW 


8.10 


0.568378 vr 2 


0.523106 vr 2 


0.62894 


TILW 


8.00 


0.561610tt 2 


0.515069 vr 2 


0.62107 


Iwasaki 


9.485 


0.420531 vr 2 


0.379397tt 2 


0.67066 


Iwasaki 


8.026 


0.420531 vr 2 


0.370379 7T 2 


0.59561 


DBW2 


12.745 


0. 153264 vr 2 


0.131857tt 2 


0.72759 


DBW2 


10.662 


0.153264 vr 2 


0.128275 vr 2 


0.67332 



Table 5: The parameters k u , k^ 1 and the average plaquette Uq for various actions and (3 
values. See eq. © for the definition of (3. The plaquette values for the TILW, Iwasaki 
and DBW2 actions are taken from [20], [2H| and [2H|, respectively. 



perturbative series in terms of the tadpole improved parameters, which must satisfy 



TI 

3" = «o*5' (58) 

st i y 

TI 

%- = «o4> i = 1 ' 2 ' 3 ' ( 59 ) 
9ti 9 

because the plaquette term in the action is a four-link term, while the other three terms 
in the action are six-link objects. These conditions are not enough to fully fix the tadpole 
improved parameters, because we have four equations with five unknowns, leaving us with 
some freedom of choice. The simplest choice is to define 

(60) 
(61) 

i = l,2,3. (62) 
However, it is important to notice that, although this definition makes the formulae for 





9 2 


9ti 


4 ' 






r TI 




°0 


= c 






c TI 


= Uq C{ 
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the tadpole improved parameters very simple, with this choice 



C* 1 = 4 1 + 8c 77 + 164 1 + 84 V 1, (63) 

which modifies many of the formulae we wrote down in earlier sections. This means, we 
have to replace every g 2 by g\ T and every c\, c 2 and C3 by cj 1 , c 2 7 and C3 7 , respectively, 
while keeping cq unchanged. The effect of introducing tadpole improved coefficients f!62|) 
is that the rescaled gluon propagator remains of the same form as we change Uq, thus 
ensuring fast convergence. 

For the perturbative calculation of Zq F we need to know the perturbative expansion 
of u 123111: 

uo = l- 9 -iffk T J, (64) 

l07T 



where k^ 1 is given by the integral 

~ w/a d 4 k 



k T J = 4vrV 



,/a (2tt)< 



kiD 11 (k,C 11 )-k 1 hD u (k,C 11 ) . (65) 



Here D n (k,C TI ) and D u (k,C TI ) are components of the gluon propagator ([T3]). with 

-3 > 



Cq = 1, C\ = C2 + C3 and C 2 = Ci — c 2 — c 3 being replaced by Cq 1 = c + 8cf 7 + 16c 2 P/ + 8c^ / 



Cf 7 = Mq Cx and Cj 7 = Mq C 2 , respectively. The numerical values of fc^ 7 are given in 
Table 03 for the various actions. For comparison, we also give the non-tadpole improved 
result (i.e. with coefficients C\ and C 2 ), which we call k u . The 'measured' values of Uq 
are also collected in Table 03 We then obtain 



= 1 - TjgF (l - "o - j ) K' ■ (66) 



Let us now introduce constants Bq(p, C) by writing the original one-loop renormalisation 
constants of Section 3 as 



2 



1 " tH" bo logM + B (p, G)\ , (67) 



167T 



where 70 is the anomalous dimension of O. We then obtain tadpole improved renormal- 
isation constants for r = 1: 



log(a/i) + S (p, C T7 ) - ( 1 - n D - - ) k. 



TI _ P U Q j 1 S^ - 

z, 



p- 4(1 -M ) t 167T 



riTI 
°0 



4 i.n 



(68) 

where B (p, C TI ) is the analogue of B (p, C), with C , Ci and C 2 being replaced by C^f 7 , 
Cf 7 and C^ 7 , respectively. In the following we shall use the abbreviation 

Bl 1 = B (p, C TI ) - ( 1 - n D - - ) k T u \ (69) 
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giving 



7 Ti 

J o 



1 — Tlr) 

P U 
p- 4(1 -M ) 



9tiCf 
16tt 2 



7o 
Co" 



log(au) + Eg 7 



(70) 



So far we have only tadpole improved the gluon propagator in our calculation of the 
perturbative factor on the r.h.s. of eq. (J4U|) . but not the fermion propagator. There is no 
need to do so for Wilson fermions. However, already for clover fermions we have seen 
that the fermion propagator ought to be improved as well. If we want the rescaled fermion 
propagator (see eq. (JEHI)) 



• u , 1 2 u o 
ip+-ap — r 

2 p — 4r(l — uq) 
to have the same form as we change Uq, we must replace p by 

TI p- 4(1 -M ) 



p 



Uq 



(71) 



(72) 



for r = 1. An alternative derivation of p TI , without expanding as a power series in a, is 
given in Appendix C. This defines 'fully tadpole improved' renormalisation constants 



vFTi 
J o 



with 



p<- nD f g 2 TI C F 

~ P -4(1 -M ) I 167T 2 

B^ = B (p TI ,C TI )- 



a 



l-n D 



ri log(au) + B F TI 



Tl 



(73) 



(74) 



Before we present numbers for Zj/ and Zq , let us make a few general remarks 
concerning the choice of gauge action. Only if g^i is a 'good' expansion parameter, can 
we expect the perturbative series to converge fast. It is generally accepted that gjjs{p) 
is a good expansion parameter for appropriate choices of u. To one-loop order we have 



2bn 



log log ■ 



A-MS 



aA 



lat 



26 log(au) + d g + Nf d 



f ■ 



(75) 



where b = (11 — 2/3A"j)/(47r) 2 , and Nf is the number of flavours. (It is appropriate to 
consider the case of general Nf here.) The ratio of A parameters is thus given by 



A 



lat 



A 



exp 



MS 



d a + d 



2bn 



(76) 



Upon inserting ()60|) and (|64p. we obtain 



2b log 



A 



log 



MS 



2b log(au) + d g + N f d f + -H- , (77) 
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Action 


Alat/ 

N f = 


Ams 
Nf = 2 


AS/ 
N f = 


Am5 

N f = 2 


Plaquette 
Symanzik 
Iwasaki 
DBW2 


0.0347 
0.184 
2.13 
45.4 


0.0172 
0.115 
1.86 
60.7 


0.380 
1.06 
5.82 
65.6 


0.262 
0.843 
5.86 
92.2 



Table 6: Ratios of A parameters. 



giving 



A lat 

Ai 



exp 



dg + Nfdj + ky/Gn 
2b. 



The coefficient d g is known for some of our actions |3Uj : 



(7S 



Action 




Plaquette 


-0.4682 


Symanzik 


-0.2361 


Iwasaki 


0.1053 


DBW2 


0.5317 



(79) 



The coefficient df is known for p = 1.4 |31j : 

= -0.01449. 



(80) 



Unfortunately, d g is not known for the Liischer-Weisz action, nor is it known for tadpole 
improved coefficients ()61)1 and (j62jl . We expect the numbers for the Liischer-Weisz action 
to be close to the result for the Symanzik action though. In Table El we have computed 
the ratios of A parameters A lat /A^j3r and A^/Ajg^ for coefficients (fT9*j) and (J%U|) and, 
to be consistent, with k^ 1 replaced by k u . We see that tadpole improvement drives Ai at 
towards A^- for the Symanzik action, giving qti ~ ffwsi so ^ na ^ 9ti appears to be a 
good expansion parameter. This is not the case for the Iwasaki action, and even less so 
for the DBW2 action. We thus may conclude that the Symanzik action, and possibly the 
Liischer-Weisz action as well, is a suitable lattice gauge field action not only from the non- 
perturbative perspective, but also from the perturbative point of view, if supplemented 
with tadpole improvement. 
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Action 


P 


Bs,p 


7 MS 


B S,P 


r?TI,MS 


tdFTI 
B S,P 


7 FT I, MS 
Z S,P 


oy iiidiiZ/iiv 


8 A^ 


— lO.O i u 


1. 


.139 










TTT W 

jl IIjVV 


r fin 




1. 


164 




1 90^ 


— ^ ^9Q 


1 1AA 


TTT W 

JL IIjVV 


R A^ 


— i i.yu i 


1. 


.168 






_q 491 

— o .^±Z 1 




TTT 

jl IIjVV O 


r ^n 


— 1 1 8R9 


1 

1 . 


I/O 


— 1 (\AA 


1 ^9^ 


— ^ ^97 


l.OOU 


TILW 


8.20 


-11.806 


1. 


.176 


-1.677 


1.337 


-3.601 


1.395 


TILW 


8.10 


-11.736 


1. 


.180 


-1.717 


1.350 


-3.683 


1.412 


TILW 


8.00 


-11.648 


1. 


.184 


-1.766 


1.366 


-3.778 


1.432 


Iwsasaki 


9.485 


-9.851 


1. 


.192 


-2.527 


1.334 


-3.851 


1.381 


Iwsasaki 


8.026 


-9.851 


1. 


.227 


-2.610 


1.481 


-4.379 


1.573 


DBW2 


12.745 


-7.270 


1. 


.354 


-4.095 


1.506 


-4.543 


1.541 


DBW2 


10.662 


-7.270 


1. 


.423 


-4.122 


1.681 


-4.666 


1.739 



Table 7: The constants B$ t p and Z<fp (an = 1) for various levels of improvement. 



Action 




By,A 


7 MS 
^V,A 


B V ,A 


7 TI,MS 
Z V,A 


tdFTI 
B V ,A 


rv FT I, MS 
Z V,A 


Symanzik 


8.45 


-13.397 


1. 


.134 










TILW 


8.60 


-11.180 


1. 


.153 


-0.099 


1.258 


-1.342 


1.290 


TILW 


8.45 


-11.100 


1. 


.156 


-0.101 


1.268 


-1.384 


1.303 


TILW3 


8.30 


-10.990 


1. 


.160 


-0.105 


1.280 


-1.431 


1.319 


TILW 


8.20 


-10.922 


1. 


.163 


-0.106 


1.289 


-1.464 


1.330 


TILW 


8.10 


-10.835 


1. 


.166 


-0.109 


1.299 


-1.410 


1.343 


TILW 


8.00 


-10.727 


1. 


.169 


-0.112 


1.311 


-1.539 


1.358 


Iwsasaki 


9.485 


-8.362 


1. 


.163 


-0.271 


1.252 


-1.196 


1.285 


Iwsasaki 


8.026 


-8.362 


1. 


.193 


-0.204 


1.356 


-1.469 


1.422 


DBW2 


12.745 


-3.538 


1. 


.172 


-0.277 


1.204 


-0.582 


1.228 


DBW2 


10.662 


-3.538 


1. 


.206 


-0.230 


1.264 


-0.607 


1.304 



Table 8: The same as Table but for By,A and Zy^(afi = 1). 
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Action 


P 


Bj' 


7 MS 

Zjrp 


B^ 


7 TI,MS 

Zjrp 


tdFTI 


7 FT I, MS 

Zjrp 


oy iiidiiZ/iiv 


8 A^ 


— 1 A 8R1 

— l^.OU 1 


1. 


.148 










TTT W 

jl IIjVV 


r fin 


— 1 9 7^0 


1. 


174 


— 1. 1UU 


1 98^ 
1 .zoo 


9 1 98 


i 


TTT W 

JL IIjVV 


R A^ 




1. 


.178 


—i ns9 


1 90^ 


—9 1 zLfi 


1 "V)A 


TTT 

jl IIjVV O 


r ^n 


— 1 9 ^71 


1 

i , 


. lOO 




1 .ouo 


—9 1 fi^ 
— z. luo 


1 .o^u 


TILW 


8.20 


-12.507 


1. 


.186 


-1.048 


1.318 


-2.177 


1.352 


TILW 


8.10 


-12.424 


1. 


.190 


-1.032 


1.328 


-2.189 


1.365 


TILW 


8.00 


-12.322 


1. 


.194 


-1.012 


1.341 


-2.202 


1.380 


Iwsasaki 


9.485 


-10.130 


1. 


.197 


-1.029 


1.279 


-1.773 


1.306 


Iwsasaki 


8.026 


-10.130 


1. 


233 


-0.806 


1.388 


-1.851 


1.442 


DBW2 


12.745 


-6.021 


1. 


.293 


-0.364 


1.210 


-0.581 


1.228 


DBW2 


10.662 


-6.021 


1. 


.351 


-0.164 


1.257 


-0.446 


1.287 



Table 9: The same as Table [7[ but for B? and Z^ IS (ap = 1). 

In Tables [?1 |H1 and |§] we compare tadpole improved and unimproved renormalisation 
constants for p = 1.4 and some selected values of /3, which are widely being used in Monte 
Carlo simulations. The DBW2 couplings (3 = 10.662 and 12.745 correspond to (3 = 0.87 
and 1.04, respectively, in the notation of The 'measured' and perturbative plaquette 
values are taken from Table El The tadpole improved renormalisation constants given 
earlier in [S] referred to unimproved coefficients Cj and unimproved parameter p. 

How good is the mean field approximation and tadpole improvement? The ultimate 
test is to compare with non-perturbative determinations of Z, which at present we can 
only do in a single case (to be discussed later). However, there is also a powerful internal 
test, which we can make with our perturbative coefficients. If the mean field Z MF of 
eqs. (J5S|) and (J57)l is a good approximation, then we should always find that the one-loop 
coefficient of Z MF is very close to the one- loop coefficient of Z^ or Zq, i.e. that 

B s ub (p, C) = B (p, C)-(l-n D - -^j k u (81) 

is small, whatever gauge action or p we use. We test this in Fig. 01 where we compare the 
original values of By with the subtracted ones defined in ([81)1 . The original coefficients 
By have large negative values, and depend strongly on the choice of gauge action and on 
the value of p. The subtracted coefficients By^ b are much closer to zero, and depend only 
weakly on gauge action or p. So we see that the mean field approximation of eqs. (J56j) 
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and (|57|) is very good at the one-loop level. If we had naively used the Wilson fermion 
result ([53)1 . this test would have failed completely, not even having the right sign. In 
Fig. El we show the subtracted coefficients in more detail. We can see that they all lie 



-5 



10 



■15 



-20 



-25 



-30 




DBW2 



Iwasaki 



T1LW 



1.2 



sub 




1.4 



1.6 



P 



Figure 4: A test of the mean field approximation. We compare the perturbative By (filled 
squares), as defined in (|H7)l. with the subtracted By h (open circles) defined in (|5Tjl . The 
subtracted coefficients are much smaller, showing that the mean field approximation is 
very good at the one-loop level. 
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1 1.2 1.4 1.6 

P 

Figure 5: The subtracted points from Fig. 0] shown in more detail. The open circles show 
results calculated with unimproved gauge actions, the open diamonds with gauge actions 
modified according to eqs. (J61j) . (J62j) and (J72I ) . All actions considered give very similar 
results. 



in a fairly narrow band, with points calculated from the original actions (solid points) 
and the tadpole improved actions (open points) in agreement. This shows how well the 
dependence of Z on gauge action and p is described by the mean field approximation, 
leaving only a small residue to be described by perturbation theory. 

As far as we can tell from the single non-perturbative result we have found for the 
TILW action at (3 — 8.45 jH]), corresponding to a lattice spacing of a ~ 0.1 fm, tadpole 
improvement drives the perturbative number closer to the non-perturbative value. But 
the discrepancy is still of the order of 10%. 

An alternative improvement scheme has been proposed in in which Z@ F is re- 
placed by the non-perturbatively computed renormalisation constant of the local vector 
current Zy P : 

ZV = ZV (¥) ■ (82) 

\ Z V/ pcrt 
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7 VI,MS 
L S,P 


yNP 


ryVI,MS 


1.478 


1.416 


1.439 



Table 10: The renormalisation constants Zgp MS (afi = 1) and Z^ I,lvia (afi = 1) for the 
TILW action at (5 = 8.45. The value of Zy% was used as input. 



rVI.MSi 



This gives 



7 np 
J v 



— ryNP 

= 



16tt 2 

Cf9ti 
16tt 2 



log(a/x) + 5q J 



5 



FT/ 



1.4 we found Z$ p ( 



ZD 



(83) 



1.416 0. We 



For the TILW action at (3 = 8.45 and p 
used this number to compute Zq in Table ITU1 This method is not applicable to operators 
with no > covariant derivatives. In that case one would have to replace the local vector 
current by an operator with an equal number (i.e. «d) of covariant derivatives. 



5 Summary 



We have computed the renormalisation constants of local quark bilinear operators for 
overlap fermions and improved gauge field actions with up to six links. The computations 
have been performed in general covariant gauge using the symbolic language Mathematica. 
This gave us complete control over the Lorentz and spin structure, the cancellation of 
infrared divergences, as well as the cancellation oil /a singularities. The price to pay is 
high. The extension to improved gauge field actions blew up the calculation tremendously. 
In some instances we had to deal with O(10 4 ) terms. Based on generalised lattice Ward 
identities we were able to show analytically that the gauge dependent part of the self- 
energy and the amputated Green functions does not depend on the choice of lattice 
fermions. 

In the limit c = 1, C\ = c 2 = c 3 = our results agree with previous calculations 
employing the plaquette action [TO! ITT] . Comparing overlap with Wilson fermions, we 
notice that the one-loop corrections have opposite sign. This is mainly caused by changes 
in the quark self-energy, and can be understood through a mean field calculation. 



We have formulated mean field (tadpole) improvement for overlap fermions. For the 
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Symanzik action the boosted coupling qti turns out to be close to gjjs-i which makes gxi 
a good expansion parameter. We thus may expect that the perturbative series converges 
rapidly. For the Iwasaki action, and in particular for the DBW2 action, this appears not 
to be true. In this case boosted perturbation theory might even worsen the situation. 

We have presented results for a variety of parameters and couplings, which cover most 
of the parameter values used in recent Monte Carlo simulations. We are happy to supply 
numbers for different choices of parameters on request. Details of our results for operators 
with covariant derivatives JSJ will be given elsewhere. 
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Appendix A 



We set a = 1 throughout this Appendix. Let us first consider the gluon propagator. We 
introduce the abbreviation 



j c n = Y i k», n = 2,4, 



The coefficients D n and D m n in eq. (jl4j) can be written in the form 



Dr. 



dr. 



D 



dr, 



m,n 



where d 



D(k 2 ) 2 ' D(k 2 ) 2 

d nt m- For D, in the case Co = 1, we find 

d = (i - d k 2 f (k 2 ) 2 - c 2 (i - d k 2 ) 2 k 2 [(k 2 ) 2 + 2 fc 4 

+ ~ C 2 2 (1 - d k 2 ) [(P) 4 + (k 2 ) 2 k 4 + 2 (A; 4 ) 2 + 2k 2 k 6 



1 
6 



24 k 10 - 7 {k 2 f k 4 + 15 k' 2 {k A f + 12 (k 2 ) 2 k 



£.6 



- 26 k A k 6 -24k' z k 



2 £.8 



54) 



15) 



?6) 



24 



and for d n and d m ^ n we obtain 
d ~- 



d 3 

do,o 

di,i 
d 2 , 2 



do,2 

do,3 
do,4 

di,2 



Dk 2 +{l-d k 2 ) 2 (P) 3 - C 2 (1 - Ci k 2 ) 7 (P) 2 P 

1 
2 

2 Pi 



+ ^C 2 



30P° + 9(P) 2 £; 



- 3 (P) 2 - 8 k 2 P + 6 k 
-25 k 4 P + 12 P((P) 2 

C 2 (1 - d P) \-2 (P) 3 + 9 P P - 6 P 



2\2 16 



4 (P) 2 P - 7 (P) 2 - 12 k k b + Uk 

1 



:C 2 



--6C 2 (l-C 1 k 2 ) (P) 2 -P 
+ 28P , 

= -6 C 2 (1 - d P) P - 2 C 2 2 U(P) 2 - 7 P 



3(P) 3 -29P P 



X5 Co 



d - (i - ci k 2 y(k 2 y + c 2 (i - Ci fc ) 



5 /c /c 2 



+ c 2 



-3(jfe 



4\2 



3 P P + 4 Jb e 



-6 C 2 (1 - d P) P - C 2 2 (P) 2 - 7 P 

4,3 : 



-14 C 2 , 



d 01 = C 2 (1 - Ci P) 2(P) 2 -3A; 



C 2 



3 P P - 5 P 



-6 C 2 (1 - Ci P) P - ^ C 2 2 [(P) 2 - 5 P 

6C 2 (1 -CiP) + 8C 2 P, 
-15 C 2 , 

6C 2 (l-CiP) + 7C 2 P. 



J7) 



(The coefficient <i mi „ should not be confused with of eq. (jUJ))- Note that D n and D m , n 
do not depend on the choice of (covariant) gauge, i.e. on £. Both D n and -D m n vanish in 
the limit G x = C 2 = 0. 



Let us now turn to the lattice Feynman rules. We omit the colour factors and the 
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gauge coupling here. We introduce the following abbreviations: 

u(p) = Lv(-p) = i sin 2 p^ + b{p) 2 , b(p) = b(-p) = ^-p 2 - p, 



w(p2,Pi) = ^^2smp2 fl smp lfl + b(p 2 )b(p 1 ) , 



(88) 

(89) 
(90) 



The Wilson quark propagator and its inverse are of the form 



S w (p) 



7 M smp^ + -p 
v 

^sin 2 p M + {^f) 



, S w \p) =i^ 7 ^sin^ + ^p 2 = X (p)+p. (91) 



The Wilson quark-quark-gluon (qqg) and quark-quark-gluon-gluon (qqgg) vertices are 
given by 



r sin 



. (P2+Pi) M , (P2+Pi) M 

-l 7 M sin h r cos 



(92) 
(93) 



with incoming and outgoing quark momenta being denoted by p 1 and p 2 , respectively. 
Note that 



X (p) +P = S w \p) = -X)P^(P'°) = -P^(P»0) • 



(94) 



The overlap quark propagator and its inverse are of the form 



2p [w(p) + 6(p)] 
For the overlap qqg vertex we obtain 



w(p) 



^(P2,Pl)= 7 / , 

Ld{p 2 )+i0(p 1 



n „ (pa , Pi ) - "^"J- V llt J (p2,Pl 



, *o(Pi) 

w(pl) 



and for the qqgg vertex we derive 



^2^(P2,Pl, fcl, ^) = V^^O^jPl) + ^22^(P2,Pl, fcl, fc 2 ) , 



(95) 



(96) 



(97) 
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with 



u(p 2 ) + w{pi) 



2(w(p 2 )+w(pi)) 



W^(p 2 ,pi,h) + W UfM (p 2 ,p 1 ,k 2 ) 



, (98) 



(99) 



where fci and fc 2 are the gluon momenta with p 1 + k± + k 2 = p 2 , and the tensor W^ v is 
defined by 



W^(p 2 ,p^k) 



1 



MP2)v 1 w ^( P2 , Pl + k) 



uj(p 2 ) + u(p! + k) u{pi + k)+ u{pi) 
x V™( Pl + k, Pl ) + V™{p 2 , Pl + fc)Xj( Pl + fc) VZ( P i + fc,Pi) 
+ V^(p2,Pi + k) VZ\P1 + fc,Pi) X (pi) 

^(P2) +W(P1 + fc) +^(pi) , Wt , i_M 
7— v— 7 ~ , w x ^0 P2 (P2,P1 + fc) XoijPl + k) 

uj(p 2 ) u(pi + k) " 



(100) 



X 



^(pi + A^i^oCPi) 



In Appendix B we need explicit expressions for the vertices contracted with the gluon 
momentum k = p 2 — p\ and k — k\ — —k 2 , respectively: 



kvr> N (P2,pi) = j2^T^pj> 

k V 2 w ' N (p, p, k, -k)k = Y^ k^V^(p, p, k, -k) k v 



(101) 



rW,N 



k v™' N {k, o) = k v™' N {o, k) , (102) 



Using the short-hand notation 

uj — uj(k) , b = b(k), X = X (k), kV[ 
we obtain 

k V™ t k V™ = uj 2 + 2pb + p 2 , k V™ Xj k V™ = 2uo 2 p + uj 2 X + p 2 X\ (103) 

and 

(104) 



k V N = ^{kv 1 w + P -Uj) . 



Differentiating X^X = uj 2 with respect to k^, and making use of V^(k, k) = —dX /dk^ 
we find 



u 2 V^(k, k) + (cos kf, + rb) sin k^ X 



(105) 
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From 



k V™ (p + k,p) = X (p) -X (p + k), 

X (p) k V™\p, p + k) + k V^(p, p + k) X Q (p + k) = u{pf - uu(p + k) 2 
X (p)X^p + k)+ X (p + k)Xl{p) =2u(p,p + k) 2 
we finally obtain 



(106) 



kV^(p,p,k,-k)k 



P 



w(p) 



us(p) 2 — u(p,p + k)' 

Ul{j)) 2 



A'o(p) 



uj(p + k) 



(107) 



Appendix B 



Gauge dependence in position space 

General argument 

We have seen in the main text that the gauge dependent terms in the Green functions 
are independent of the gauge action and of the fermion action: they do not depend 
on the coefficients c\,C2,cs and on p, and are the same for overlap fermions and clover 
fermions [2l)j . The numbers are, however, different on the lattice and in the continuum, so 
this is not simply due to an ultraviolet convergent integral, where the lattice regularisation 
plays no role. The independence from the gauge action is fairly easy to explain, but the 
independence of the fermion action is a remarkable observation, and is something which 
we ought to understand. 

In this subsection we give a brief explanation of this phenomenon in configuration 
space, while in the next subsection we give a more formal argument using the Feynman 
rules. 

At the one-loop level we can split the functional integral over the gauge fields into 
an integral over the (physical) transverse gauge fields and the (unphysical) longitudinal 
gauge fields. The integral over the longitudinal fields is a simple Gaussian integral, given 
by the gauge fixing term added to the action. The contribution of the longitudinal modes 
is a universal factor multiplying the Landau gauge Green function 



G(x,g 2 ,0 = G(x,g 2 ,{i = l) 



1-(1- i)g 2 C F 



d A k 1 — cos kx 



(27T) 4 (#2)2 



+ ■ ■ 



(108) 
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where G can be any quark Green function, e.g. the quark propagator or any of the three- 
point functions we consider in this paper. See Appendix A of Ref. [31], where this is 
proved to all orders for non-compact QED. It is not clear whether this result will hold to 
all orders for a non-Abelian group, but it is certainly true at the one-loop level. 

Equation (jl()8j) gives us an exact expression for the 0(g 2 ) gauge dependent term in 
configuration space 

G-gauge/ ) = 2 C £ gtree, ) f 1 ~ COS ^ ^ 

K J 1 ' J (2vr) 4 (jfe2)2 1 ; 

When we Fourier transform this expression to find the gauge term in momentum space, 
the integral is dominated by x 2 ~ 1/p 2 - Since we are interested in a 2 p 2 <C 1, we have 
to consider what happens to eq. (jl(J9p when x 2 ^> a 2 . Although the tree-level Green 
functions do depend on the fermion action when x 2 ~ a 2 , at large distances all fermion 
actions agree: they simply tend to the continuum value G^ t (x). Since this is the region 
that dominates the Fourier transform, we see immediately that the gauge term will be 
the same for all fermion actions (or to be more precise, that differences in the gauge term 
due to differing choices of fermion are suppressed by powers of a). This completes our 
argument for the universality of the gauge dependent terms. 

Note that although the gauge dependent terms are independent of both gauge and 
fermion action, they would change if we used a different lattice discretisation of the gauge 
fixing term (d^A^) 2 , as advocated in 

Calculating the gauge dependent contributions 

Equation (jl(J9J) not only shows that the gauge dependent terms are universal, it also 
allows us to calculate them. At large distances (x 2 ^> a 2 ) the integral in eq. ()109|) takes 
the value 

f d 4 k 1 — coskx 1 x 2 \ , -i \ 

JWVHW * l6^( F » + ^ + '° g 4^)- < U0 > 

In position space the massless quark propagator S^ t (x) and the three-point function 

-ftrcc . 



^ccmt( a; ) are 



,1 nr, ,., 1 



Otree/ \ 

°contW 27r2 ^ 2)2 , 

-<tree / \ _ r F ' 



GSS(x)=7pT i '7a 



in 2 



2- 



*En*Err Of) 



J p tXj (7 w p(J 



X 2 ) 2 X 2 



Let us first calculate the gauge term in the propagator. We have argued that this will 

be 



2 ' i 



16tt 



(112) 
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where T is the Fourier transform. Some details for performing T are given at the end of 
this subsection. From 

^raoo] =h ( n3 ) 



i 

and 

T 



J_[l_2 TB -l g(aV)] (114) 



we get 



S( P , g 2 , = S{p, g 2 , e = 0) + ^ £ [F - 1e + 1 - log(aV)] , (115) 

which agrees with eqs. ([3*T|). (|52*|) in the body of the paper. 

Now we calculate the gauge dependent parts of the non- amputated three-point func- 
tions: 



(^o + IE + log ^) ] = ± [F - 1E + 2 - log(aV)] + i , 

(116) 



which leads to 



G'W.O = G'W.O) + ^ [F - lE + 2 - tog(aV)] + , 

(117) 

where Y x is defined in (|S5|l . When we amputate this using the gauge dependent quark 
propagator ()115|) we find 

A*(p, g\ = A x (p, g\ 0) + ^£ jr* [-F + lE + log(oV)] - ^ #7^7, * 

, (118) 

valid for general T x . The final term simplifies when we put in specific choices for T x , 
giving 

A S ' P (P, 9 2 , = A*' P (P, <? 2 , 0) + [-F + 7^ - 2 + log(aV)] {1, 75} , 

A^(P, <r\ = A^(p, 0) + ^ {7m [-Fo + 1e - 1 + log(oV)] + 2^} {1, 75 } , 
A^(p, <? 2 , = A^(p, <7 2 , 0) + [-Fo + 7£ + log(oV)] <V75 , (H9) 
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which reproduces all the gauge dependent terms in eq. (J42|). Equation (jT09j) still holds if 
a quark mass is included, so it should also reproduce the mass dependence of the gauge 
terms previously calculated in [26 . 



The integrals needed to perform the Fourier transform into momentum space are easily 
found by introducing an additional integration over a, to make the x integrations into 
Gaussians. As an example, let us transform S^ t {x){x 2 / Aa 2 ) 6 into momentum space: 

1 r , r M / r 2 \ S 1 f , r°° do o 1 ^ 6 

7' 



2tt 



^( 4a2 '" S f«fa„-« e ^ (120) 



4 r(2-5) 



p 2 r(2 - 5) 

Expanding both sides to first order in 5, we get the result ([114)1 . The integrals needed 
for the Fourier transform of the three-point functions can be obtained in the same way as 
outlined above. We obtain 

/ ^ (l + * log ^) = ^ {1 + ^ [-27* - log(aV)] } (121) 



and 



/ (l + <nog£^) = 2^ {1 + 5 [1 - 2 7 , - log(aV)] } 

(p 



47r 2 ^{l + < 5[2-2 7B -log(aV)]} 



:i22i 



Gauge Dependence in Momentum Space 

Here we calculate explicitly the gauge dependent one-loop contributions to the quark 
self-energy and the amputated Green functions using the Feynman rules. We now set 
a = 1 again. 

The generalised lattice Ward identity 

As can be easily checked, the generalised lattice Ward identity for Wilson and overlap 
fermions to lowest order is of the form 

S W ] N (P2) - S w ] N ( Pl ) = (P^2),vZ' N (P2,Pi) = (p7 r p~2)Vl v > N ( P2 ,p 1 ) , (123) 
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where inverse propagators and vertices are defined in Appendix A. ^From this identity it 
can be seen that 



t = -k V x ' (p,p + k) S w ,n{p + k) + Sw N (p)S w ,n{p + k) 
= ~S w ,n{p + k)k V™ ,N (p + k,p) + S w , N (p + k)S^ N {p) 
which leads in the limit p — > to 

1 = -k V^ N (0, k) S W)N {k) = -S WjN (k) k V^ N {k, 0) . 



(124) 



(125) 



In addition we have the ordinary lattice Ward identity to lowest order (with zero gauge 
boson momentum) 



d 



SwAv) = ~V^ N (p,p) . 



(126) 



Differentiating (j 124)1 with respect to p M , and taking the limit p — > 0, we get 

^-s w>N (p + k) k v™' N {k, o) = s w , N (k)^- k v™> N ( P + k, P ) + iS w>N (k) 7At . 

opf, P =o dpf, P =o 



Using Sw N Sw,n = 1 and ()126|) . ()127J) . another useful form can be derived: 





kV™> N {p + k,p) =Vr/\k,k) + i 7ll . 



W.N i 



p=0 



(127) 
(128) 



Gauge dependent one-loop corrections 

We consider the gauge dependent part of the gluon propagator in a general covariant 
gauge, 

D^ c (k) = ! ' J " 



(129) 



and use the short-handed notation for the D-dimensional one-loop integration variable 

f _ f dPk 
J k -J (2^- 



(130) 



For finite integrals we replace the dimension D by four. The basic divergent lattice integral 
(in dimensional regularisation) is 



+ F - log4vr - log(a 2 // 2 ^ 



! k (k 2 ) 2 16tt 2 LD - 4 
Furthermore, we use the finite lattice integral 



— = z . 

kk 2 



(131) 



(132) 
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First we demonstrate that the gauge dependent contribution of one-loop corrections to 
local quark bilinears does not depend on the particular representation (Wilson or overlap) 
of the lattice fermions. Since the local operators do not contain gluon operators, we have 
to consider only the vertex correction to the amputated Green function: 

I x (p) =g 2 C F J2 [ V^ N { P , P + k)S w ^ N { P + k)T x S w ^ N { P + k)V^ N {p + Kk)Df:nk). 

(133) 

The corresponding correction is ultraviolet (UV) logarithmically divergent. Using the 
technique for analytic handling of the divergences, we have to consider 

I x ^(p) = I x (0) + I Xc ™\ p ) . (134) 

With the propagator (|129jl . and using eq. (|125j) . we get at zero momentum the expected 
result 

rX/n\ _ „2^y / 1 Lt/-W,N, n , „\ n / K\ tX c tu\ui/W,N, 



i x (o) = - g 'c F / -^—kvr^ (^k) s w ^k)v x s w ^k)kvr^ {kq) 

'k (k 2 ) 2 

(135) 

-g 2 C F Y x 1 



(k 2 ) 2 

The pole cancels with that of the continuum contribution, I Xcont ( P ), while the one- loop 
correction is independent of the form of the lattice propagator as result of the Ward 
identity. 

For operators with derivatives one has to consider the Taylor expansion of the lattice 
integrals up to the corresponding order of the UV divergence. In that case additionally 
the 0(g) and 0(g 2 ) contributions of the operators have to be considered. It is not difficult 
to perform a similar proof. 

The sunrise diagram (i.e. the left diagram of Fig. [2J) for the quark self-energy is of the 
form 

/— (p) = g 2 C F [ V^ N (P,P + k)S w , N ( P + k)VZ' N (p + k,p)D^°(k) . (136) 

Since that one- loop integral is UV linearly divergent, we have to calculate (using dimen- 
sional regularisation) 

jsunrisc^ = pmnrise/g) + ^ # jsunrisc t p \ + Jcont/ p N _ /jgyx 

d Pa P =o 
For the finite tadpole contribution (D=4) we have 

/tadpole (p) = g 2 Cp J- f V^ v N ( P ,p, k, ~k)D^(k) . (138) 
ixv Jk 
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The gauge dependent quark self-energy contribution I lat (p) is then given as sum of 
and ()138|) . To show the independence on the lattice fermion representation, we transform 
the difference 7 lat (p) - J cont (p), using eqs. I)125|) - l)128|) . to the form 



I lat (p)-I cont (p)=g 2 C F 



where 



J^ N (t) 



2\2 



{J W,N + jW,N + jW, N) 



(139) 



kV™' N (k, 0) 



(A' 



2\2 



k (fc 



2\2 



kV^' N (p,p,k,-k)k. 



k (k 



2\2 



(140) 
(141) 

(142) 



It remains to show that the sum of integrals J^' N + J^ ,N + J^' N vanishes. For Wilson 
fermions we immediately obtain 



jr(i) = ^o, jy> 



(143) 



Therefore the sum is zero. 

Integrals for overlap fermions 

In the following we use the abbreviations (|1U2|) . Using (jl(J4j) . and taking into account 
the symmetry in the integration, we get immediately 



jf (11) = p 



1 LU + b 



k (k 



2)2 UJ 



(144) 



To calculate we use the symmetric part of V{t(k, k) in k and get 



N i 



J k 



1 1 



(k 



2\2 o;3 



uj 2 {uj — pcosk^) + p sin 2 k^cosk^ + rb) 



(145) 



To calculate J% , we split it into two pieces, J* = + Jg, according to eq. flH7|). After 
symmetrisation and taking ^ V£[ (p, p) in the zero lattice spacing limit, we find for the 
k independent part 



k (k 



2\2 



kV»(p,p)k = ij-[l-4-)Z . 



(146) 
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For the cancellation to be shown, we write this result in the form 



J3M = iXX 7 ^ 



2\2 



1 — cos kfj, — 



b+_p 
P 



(147) 



The k dependent vertex part arising from ()99|) contributes both to the unit matrix 1 and 
to p'. After integrating over k we have to take the zero spacing limit. This is achieved by 
performing a Taylor expansion around p = up to linear terms in p: 



2\2 







k (k 

J£(l) + Jg( 



k v£ (o, 0, k, -k) k + ^p^—kvg (p, P , k, —k) k 



p=0 



Using the form (|107|) we get 

1 ijj + b 



-p 



j. 



32 ( 



k (k 



2)2 LU 



k (k 



2\2 



b + p (u — p) cos k^ p sin 2 k^ (cos k^ + rb) 



P 



0J 



(148) 



(149) 



We see that J32W cancels the contribution from the sunrise diagram ()144j) . and 



0. 



(150) 



6 Appendix C 

In this appendix we give a more complete derivation of the mean field overlap fermion 
propagator and tadpole improved p, without having to expand in powers of a. 

To calculate the mean field overlap fermion propagator, we start from the mean field 
value of the Wilson-Dirac operator in momentum space, 

Dyrjug, r) = - f up ^ i 7m sin a P» + 4r - u r y~] cos ap^ J . (151) 

V A* A» / 

This can be written in terms of the tree-level Dirac operator 

4r 

D™ F (u , r) = u D^ c (r) + -(1 - u ) . (152) 
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Let us now calculate the overlap operator from D^ F - The first step is to find the X 
operator 



X^(u 0: r,p) = D^(u 0: r)- 



a 



= u 



D w( r ) 



p - 4r(l - m ) 



au 



where 



= u X^(r,p TI ), 

TI _ p - 4r(l - up) 
u 



We are now ready to calculate the mean field overlap operator 



D^(u ,r,p) = ^ 

(X 



X (u , r , p) 



P 
a 



X MF \u ,r,p)X MF {u ,r,p) 

X trec (r,p r/ ) 
sjx^ c \r, p TI )X^ c {r, p TI )_ 
9 DT(r,p TI ) = Zf F D%™(r,p T >) 



1 + 



(153) 



(154) 



(155) 



In other words, the mean field overlap operator is proportional to the tree-level overlap 
operator, calculated with the same r but with a tadpole improved p. The constant of 
proportionality is 

7 mf P . P u o 



p TI p _ 4r ( 1 _ Mo ) 



(156) 
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